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We study the quantum behaviour of chaotic billiards which 
exhibit classically diffusive behaviour. In particular we con- 
sider the stadium billiard and discuss how the interplay be- 
tween quantum localization and the rich structure of the clas- 
sical phase space influences the quantum dynamics. The anal- 
ysis of this model leads to new insight in the understanding 
of quantum properties of classically chaotic systems. 

PACS number: 05.45.+b 



I. INTRODUCTION 

When in a cold winter morning at the start of 1976 one 
of us (GC) entered the Institute of Nuclear Physics of the 
Siberian Division of the USSR Academy of Sciences in 
Akademgorodok, the main purpose was to discuss with 
the author Boris V. Chirikov, a startling paper which 
appeared as a CERN preprint at the end of the sixties 
jl| . Except to few, the paper remained almost unnoticed 
and practically none had the perception that this paper 
was a cornerstone of the new building of nonlinear science 
which, only few years later, expanded in an impressive 
way. Among other things, the so-called standard map 
or Chirikov map was discussed as a model of classical 
dynamical chaos. It was in that same occasion that we 
quantized the map thus obtaining the quantum kicked 
rotor which became later a basic model for the study of 
quantum chaos [§,1). 

The importance of the Chirikov map lies in the fact 
that it gives an approximate description of a wide class of 
dynamical systems. As a matter of fact it contains most 
of the complex features of dynamical motion. For these 
reasons its properties are continuously being studied and 
it is used in university courses to teach the emergence of 
random motion in deterministic systems. 

In the quantum world the standard map played a fun- 
damental role in bridging different fields of research: solid 
state physics, atomic physics, Random Matrix Theory. In 
particular it allowed to discover the phenomenon of quan- 
tum dynamical localization namely the fact that quan- 
tum interference effects suppress the classical diffusive 
process and lead to localization of the quantum excita- 
tion. This phenomenon had then been observed in several 
laboratory experiments in atomic physics and, more 



recently, a physical microscopic realization of the kicked 
rotator had been possible J8|. 

All the above experiments refer to systems under exter- 
nal perturbations. It was however not clear how quantum 
localization could take place in conservative systems. In- 
deed in this case an Eigenfucntion (EF) is always a super- 
position of a practically a finite number of unperturbed 
EFs and therefore it is not even clear what localization 
should mean here. 

In Refs. [p|jic|l, by using a well known model of quan- 
tum conservative systems, the so-called WBRM model, 
it was shown that EFs can be localized inside the energy 
shell. More precisely, an ergodicity parameter A was in- 
troduced which acts as a scaling parameter. For A < 1 
the EFs are localized inside the energy shell 6E, while 
for A > 1 the EFs are extended or ergodic. The width 
8E is the width ot the Local Density of States (LDOS) 
or strength function. Therefore for A > 1 an EF actu- 
ally excites the maximum number of unperturbed states. 
As it is quite evident this problem calls in question the 
notion of quantum ergodicity JTT1 . 

It was then shown that quantum localization actually 
takes place in a more physical model of conservative sys- 
tems: the billiard in a stadium [Q. The same effect 
was successively confirmed in a different shape, classi- 
cally chaotic, billiard 13 14|]. 

The question remains open however, whether or not 
the mechanism of localization is the same as in WBRM. 
In particular it is necessary to understand how the struc- 
ture of EF changes as one moves from the perturbative 
regime to the region of complete ergodicity. 

In this paper, on the basis of the results of M] , we dis- 
cuss in detail the above problem. We show that contrary 
to WBRM, localization here manifests via the sparsity of 
the structure of the EFs and we surmise that this should 
be a general feature of localization in conservative sys- 
tems. We also show that the presence of cantori in the 
classical motion plays an important role in the quantum 
motion. The "quantum cantori" should be observable in 
real laboratory experiments. 



II. THE CLASSICAL STADIUM BILLIARD 

In this paper we mainly consider a stadium billiard, 
that is the motion of a free point particle of unit mass 



and velocity v (energy E — v 2 /2) bouncing elastically 
inside a stadium-shaped well: two semicircles of radius 
1 connected by two straight line segments of length 2e. 
The classical motion, for arbitrary small e, is ergodic, 
mixing and exponentially unstable with Lyapunov expo- 
nent A ~ e 1 / 2 . It can be approximated (up to 0(e)) with 
the discontinuous stadium-map Jl2l 



L n+ i = L n - 2esin6»„sgn(cos6»„) v /l - £ 2 , 
9 n+ i = 8 n + tt — 2 arcsinL n+ i (mod 2ir) 



(1) 



where L = lj \J2E is the rescaled angular momentum, 
/ = r A v, and 9 is the polar angle (identical to the arc- 
length for small e). Note that, due to symmetry, the 
stadium map ([!]) has a period tt in the angle 8. One 
should note also that the stadium map is topologically 
similar to the sawtooth map 



Pn+l = Pn + K(x - |) 

x n+ i = x n +p n +i (mod 1) 



(2) 



The behaviour of the sawtooth map (2) has been studied 
analytically and numerically in great detail [^^7). For 
the stadium billiard, due to the classical chaotic motion, 
the angular momentum I undergoes, for small e, a process 
of diffusive behaviour 



((In ~ h) 



D(l )n 



(3) 



with diffusion rate D — D(l) which may depend on local 
value of angul ar mom entum I. Rigorous results on the 
have shown that the diffusion rate 



V, 



sawtooth-map |16 

goes as D cx K^UK < 1 or D cx K 2 if K > 1. The 
same dependence D cx e 5 / 2 on the small control param- 
eter e has been numerically found in |l^ | for the billiard 
and for the map (0). The factor y/l - L 2 n in the local 
shift of angular momentum leads to a dependence of dif- 
fusion rate on the local value of angular momentum. As 
a result, it is easily seen that, for the map ([!]): 



D(l) ^ 1 e 5/2 (2E-F), 



(4) 



where the prefactor 7 « 2 has been numerically com- 
puted. In figure 1 we show the weak dependence of 7 
on the local angular momentum I for various values of 
the parameter e in the stadium billiard. Note that eq. 
(||) gives the diffusion rate in the discrete time n. The 
physical time T between collisions depends on the lo- 
cal angular momentum, namely T = \/2E — l 2 /E. Since 
the stadium billiard is mixing, the classical dynamics will 
lead to the microcanonical equilibrium distribution. This 
implies that the equilibrium distribution of angular mo- 
menta is uniform in the discrete time 

PdiscrCO = 1/V8E 

(invariant measure of the approximate Poincare map ([j])) 
while, in the continuous time, it obeys the semicircle dis- 
tribution 



p(l) = ^-V2E-PocT(l)p d 
nE 



(5) 



(invariant measure of the smooth dynamics in the full 
phase space). 

The power 5/2 in eq. ([|) (instead of 2 as found for 
smooth deformations of the circle [^3|), is due to the ex- 
istence of cantori in the classical billiard motion which 
are also typical of discontinuous maps like (0J^) The can- 
tori form strong obstacles to phase space transport and 
therefore they reduce the diffusion rate (even if the dif- 
fusion remains normal). 

It is important to evaluate the flux T of the phase 
space area which is transported through a cantorus in 
one iteration of the map. For the sawtooth map there 
are no dominant cantori in phase space and the flux T is 
independent of the winding number of the resonance and 
is given by T = K 2 j%\[D where D = K 2 + AK g§ . 

Therefore, for small values of the control parameter K 
the flux goes as fsT 3 / 2 . We can now use this result to 
estimate the flux through cantori for the stadium map 
for small e. We obtain 



T w (2E) 



l/2 £ 3/2 



(6) 



which incorporates the correct energy scaling of phase 
space area. 

The size of a cantorus C, in the rescaled angular mo- 
mentum variable L, that is pc — max^ L — mine L, av- 
eraged over all the resonances, can be estimated from 
the exact results on sawtooth map . After averaging 
over all the resonances one finds that the average size is 
proportional to the parameter e, 



P = ce, 



(7) 



where c is some numerical constant. From our numerical 
computations on the stadium, c turns out to be in the 
range c = 10 for cps=0.05, and c = 15 for eps =0.005. 
The fact that c slowly increases with decreasing eps is 
due to existence of the cantorus along the separatrix of 
2:1 resonance (around L = 0) which has a larger size 
and scales as p(2, 1) w y/e. In figure 2a we illustrate the 
structure of cantori by showing a phase space portrait of a 
typical orbit (e = 0.003) around the largest 2:1 resonance. 



III. THE QUANTUM DYNAMICS 

In this section we consider the structure of quantum 
eigenfunctions for the stadium billiard v I'„(r) which are 
solutions of the Helmholz equation 

(A¥„ + = 0. 

with consecutive eigenenergies E n = k 2 /2. We consider 
only odd-odd states with Dirichlet boundary conditions 
on the quarter-stadium. It is interesting to consider the 
eigenstates of the stadium billiard expressed in terms 



of eigenstates of the nearest integrable billiard, namely 
the circle billiard. The eigenfunctions of the unit quarter- 
circle billiard are 

$ sm (r) = J2 S (k° sm r) sm(2s4>), 

where k® m are zeros of even-order Bessel functions. One 
may expand an eigenstate of the stadium billiard (for 
small [j] e) in terms of eigenstates of a quarter-circle bil- 
liard 



sm 



(8) 



The probability of having a value I — 2s of angular mo- 
mentum (only even values of angular momentum are al- 
lowed due to symmetry) is 



p n (l = 2s) 



\ C sm\ 



(9) 



If quantum eigenstates were ergodic one would expect 
to recover the classical microcanonical distribution (||) of 
angular momentum namely one would expect the quan- 
tum distribution 



Pn {l = 2a) w p(l) oc y/2E n - P. 



(10) 



However, since the classical motion is approximated 
by the map (1), then in analogy with the quantization 
of the standard map, one would expect the phenomenon 
of quantum localization to take place with localization 
length £ proportional to the classical diffusion rate 



D. 



(11) 



The localization length £, which is defined more precisely 
below (eq. (16), section 6), measures the (average) num- 
ber of excited angular momentum eigenfunctions, that is 
the average width of the probability distribution p n (l). 

In a bound conservative system like the stadium how- 
ever, there are two main peculiarities which influence the 
process of quantum localization. First, the range of varia- 
tion of angular momentum / is finite, —\[2~E < I < \J2E. 
As discussed in [jl2|, this implies the existence of an up- 
per bound for the localization length which is defined 
by the condition € max = ^max = \2E. This leads, to- 
gether with ([ll]), to the derealization or ergodicity bor- 
der e 5 N ~ 1, where N w E/8 = fc 2 /16 is the sequen- 
tial quantum number in the stadium billiard. Therefore, 
in order to observe 'dynamical' localization, this border 
must be above the perturbative border. The perturbative 
border is given simply by the condition ey2E ~ 1 , which 



If e is not really small, in order to be able to express the 
eigenfunctions of the stadium in terms of eigenfunctions of 
the circle, one should consider the smallest enscribing circle 
billiard with radius 1 + e. 



states that the length of the straight line segment of the 
stadium should be at least one de Broglie wavelength in 
order to be visible to quantum mechanics e ~ (2EY~ 1 ' 2 . 
The same border can be obtained from the map (pj) via 
the condition that in one iteration the change in angu- 
lar momentum is at least one quantum Al ~ e\/2E ~ 1. 
Therefore, for small values of the parameter e the two 
borders are well separated 



N v = —e 
p 16 



« N > = ST" 



and is possible to observe localization as actually done in 
S 

The second peculiarity is due to the presence of can- 
tori. Already long ago it has been surmised by McKay 
and Meiss [jljj that cantori may act as perfect barriers for 
the quantum motion if the flux through cantori, for each 
iteration of the quantum map, is less than one Planck's 
cell. However, to our knowledge, this effect has never 
been observed neither numerically nor experimentally 
and therefore it is not really known whether the above 
condition actually plays a role in quantum mechanics. 
For our present case of the stadium map this border can 
be estimated from the flux (|) T ~ x := e 3/2 V2E ~ 1 
(note that we have h = 1), and gives 



N c = — e-\ 
16 



(12) 



It is interesting to observe that, for sufficiently small e, 
this border is well separated from the other two borders 
N p <C N c <C N e and should be (numerically or experi- 
mentally) observable. Hence, for small values of the con- 
trol parameter e, we expect to see four different regimes in 
the quantum dynamics of the stadium biliard: (1) per- 
turbative regime for N < N p , (2) pseudo-integrable or 
cantori regime where eigenstates are expected to be lo- 
calized on classical cantori, for N p < N < N c , (3) dy- 
namical localization for N c < N < N e , and finally (4) for 
sufficiently large energy, N > N e , we should enter the 
regime of quantum ergodicity where (10) holds. 



IV. CANTORI AND QUANTUM MECHANICS 

Here we would like to illustrate explicitly the effect of 
cantori on quantum eigenfunctions. In the regime where 
x = e 3 / 2 fc < 1 it is natural to expect that cantori will in- 
fluence the localization process. Indeed, as shown above, 
in such situation the cantori act as perfect barriers, and 
the quantum system looks as if classically integrable [^0) . 
It is therefore expected that the localization length of 
eigenstates must be of the order of the size of cantori. 
On the other hand, since for small e the cantori border 
can be much higher than the perturbative border, we 
may have here a nice possibility to study the effect of 
cantori in quantum mechanics. The average size of can- 
tori has been found analytically and numerically to be 
(see section II, eq. (^)) 



4 = ceZ max = ceV2E (13) 

with c a numerical constant. In figure 2b we show an 
eigenstate in angular momentum basis in the cantori re- 
gion near the 2:1 resonance. In figure 3 (and also figure 
4) our numerical data show that at x = 1 the behaviour 
of localization length changes: for x < 1 namely below 
the cantori border the localization length agrees with the 
theoretical estimate (|~3|). This provides the first numer- 
ical evidence that the idea of flux quantization through 
cantori introduced in |l9[ is indeed correct. Even more 
convincing is the inspection of individual quantum eigen- 
states in the Husimi phase space representation (see fig- 
ure 5) which show very clearly, that in the cantori re- 
gion, the relative phase space area occupied by the lo- 
calized eigenstates does not increase with increasing e or 
with increasing k = \J2E. We checked this even quan- 
titatively, by computing the phase space area of Husimi 
phase space distributions via information entropy or in- 
verse participation ratio p2] | . The width of Husimi func- 
tions on quantized cantori in the perpendicular direction 
is given only by the width of the wavepacket which is used 
for the computation of Husimi functions (see section VI) . 

V. INTERPLAY BETWEEN CANTORI AND 
DYNAMICAL LOCALIZATION 

Above the quantum cantori border x — 1 the flux 
through the turnstiles becomes larger than a Planck's 
cell, the cantori do not act any more as barriers for quan- 
tum dynamics and the quantum motion starts to follow 
the classical diffusive behaviour up to the quantum relax- 
ation time tft (break time), which is proportional to the 
density of operative eigenstates namely of those states 
which enter the initial condition and therefore actually 
control the quantum dynamics. For t > tn instead, the 
quantum dynamics enter an oscillatory regime around 
the stationary localized state. The density of operative 
eigenstates is by a factor ^/Z max smaller than the total 
density of states 

t R = (l/l max )dN/dE 

and therefore the relaxation time is less than the Heisen- 
berg time which is given by the level density. The angular 
momentum width i of the localized state is then given by 
I 2 w Dtn/T, where T w E~ x l 2 is the average time be- 
tween bounces. This leads to the simple expression for 
the average (scaled) localization length 

£/l max « D/k = ae 5 ' 2 k (14) 

where a is a constant to be determined numerically. We 
need however to take into account the fact that for x < 1 
the width of eigenfunctions is determined by classical 
cantori and not by dynamical localization. Since we mea- 
sure the localization length in angular momentum space, 



then for x > 1, we need to add the average size of cantori 
to the angular momentum spread due to dynamical lo- 
calization. Therefore, above the cantori border x > 1, we 
expect the following expression for the numerical (scaled) 
localization length a n — £ n // max , 

cr n = P + (1 - p)ae(x - 1) (15) 

which takes into account the fact that we need to rescale 
the total size of angular momentum space, and that for 
x = 1, cr„ = p which is the average size of cantori as 
determined numerically. In figure 3 (see also figure 4) 
it is seen that expression ( p"5| ) is in excellent agreement 
with our numerical data up to N ~ 10 and for different 
values of e. The obtained value of the numerical constant 
is a = 1.7. 

Notice that only for verysmall e the constant contri- 
bution due to cantori in ( |15| ) will be negligible and lo- 
calization length would be simply proportional to D. as 
given by eq. (pT|). 

Notice also that due to the discontinuity of the stadium 
map, the quantum eigenfunctions are not exponentially 
localized, like for the standard map or for smooth chaotic 
billiards |ll| , but they have power-law tails 

p(l)« \l-lo\~ 4 

This is consistent with results based on the quantization 
of the stadium map pl[ , and with rigorous results on 
band random matrices with increasing band size p3| . In 
figure 6 we show two (averaged) localized states in an- 
gular momentum basis in which the power law tails are 
clearly seen. 

The pseudo-integrable (cantori) regime may also be de- 
tected by inspecting the energy level statistics, e.g. the 
commonly studied nearest neighbour level spacing distri- 
bution P(S), the delta3 statistics etc. As an example in 
figure 7 we show that P(S) is nearly Poissonian in can- 
tori regime while it is intermediate between Poisson and 
Wigner (GOE) in the regime of dynamical localization. 
We would like to stress that the above deviations from 
GOE predictions have no relation with periodic orbit the- 
ory and bouncing ball orbits. 

VI. LOCAL DENSITY OF STATES 

The Hamiltonian H of the stadium billiard may be 
written as a small deformation of the integrable cir- 
cle billiard Ho, namely H = Hq + eV. The exact 
stadium-eigenstates may be expanded in terms of unper- 
turbed quarter-cirlce eigenstates $ sm (r) where k° sm are 
the eigenvalues of the integrable quarter-circle - the ze- 
ros of the even-order Bessel functions (@). Or, vice versa, 
we can express eigenstates of integrable quarter-circle in 
terms of exact eigenstates of the perturbed billiard by 
the same (orthogonal) matrix of coefficients c™ m , 

J I 



It is important to note that unlike for Wigner-band- 
random matrices [Tc| l the matrix of coefficients c™ m , or- 
dered with increasing wavenumber (energy), has been 
found to have a symmetric appearance (see figure 8a for 
an example). The structure of rows (expansions of ex- 
act states in terms of circle states) is very similar to the 
structure of columns (expansions of circle states in terms 
of exact states). The effective bandwidth b of the ma- 
trix c" m determines the width of the energy shell SE = b 
while the effective number of nonzero entries in each row 
(fixed n) is proportional to the localization length £ of 
that state fy n . Such a symmetry between perturbed and 
unperturbed states seems to be generic for conservative 
Hamiltonian systems. Indeed, a similar structure of the 
matrix of coefficients has been found for the chaotic rough 
billiards (with shapes of wiggled circles) introduced in 
Jl3| as well. However, the bandwidth b for the stadium 
billiard has been found to take always its maximal value, 
that is 



independent of the parameter e, except in the pertur- 
bative regime N < N p where it becomes smaller. This 
should be contrasted with chaotic rough billiards |13| for 
which, below the so called Wigner-ergodicity border, the 



bandwidth b has been found to decrease as b ~ e N. 

Below the ergodicity border, N < N e , i < Z max , the 
matrix c^ ns is sparse in both horizontal and vertical di- 
rections (see figure 8a). (Note that the same is true also 
for the rough billiard) As a consequence of that, the av- 
eraged local density of states (LDOS), 

w(k) = (J2\CsJ 2 6(k + k sm -k n )) sm 



(averaged columns shifted to the same center) is nearly 
the same as average EF in circular basis 



W(k) = <]T KJ 2 S(k + k n - (n\^2H Q \n))) n 

sm 

(average rows shifted to the same mean). See figure 8b. 
Both distributions, W(k) and w(k), agree quite well with 
the theoretical distribution W e (k) — (sin(k) / k) 2 / tt which 
has been proposed |l3|] for LDOS of nearly circular bil- 
liards in the regime of quantum (Shnirelman) ergodicity. 

Such structure is found in the regime of dynamical lo- 
calization as well as in the regime of cantori localization 
where the quantum system behaves as if classically in- 
tegrable. We would like to stress that for the case of 
WBRM there is instead a strong asymmetry in the sense 
that the EF's are solid and narrow, namely they are lo- 
calized inside the energy shell, while the expansion of a 
basis state in terms of exact EF's shows a sparse struc- 
ture [El. 



VII. NUMERICAL METHODS 



A. Computation of eigenstates of quantum stadium 
billiard 



In order to check the above theoretical predictions we 
had to compute numerically the quantum cigenfunctions 
^n(r) = of chaotic billiards which are solutions 

of the Schrodinger equation (A + fc J 2 J'i , n = (where 
Ti = 1) satisfying Dirichlet boundary conditions along 
the boundary. We had to compute eigenstates with 
extremely large sequential quantum numbers N of the 
desymmetrized stadium billiard (and also of the rough 
billiards) up to N w 10 7 . 

We have used a recently proposed scaling method Q , 
but we expanded quantum eigenfunctions \&fc(r) in terms 
of more suitable circular waves (here we consider only 
odd-odd states) 



M 

E 



a s J2 S {kr) sin(2s#) 



instead of the originally proposed plane waves. The 
eigenvalue k = \/2E is determined by minimizing a cer- 
tain positive bilinear form along the boundary of the 
billiard |p4| . This can be done by solving a general- 
ized eigenvalue problem of dimension sa k which yields 
~ 10%/c accurate eigenvalues k n and eigenvectors a™ si- 
multaneously. Therefore, with this method, the com- 
puter workload required per energy level is by a factor 
k or even lOfc (namely from 10 3 to 10 4 for k = 10 3 ) 
smaller than with more traditional methods, like the orig- 
inal Heller's plane wave decomposition or boundary in- 
tegral method! 



B. Definition and computation of localization length 

Note that the magnitude of coefficients a™ is propor- 
tional to the probability (|^) of having angular momentum 
I = 2s in a quantum state \t n , 



p n (l = 2s) = \a s \ 2 / drr\J 2s (kr)\ 2 cx \a s \ 2 ^ 1 - I 2 / k 2 . 



The number of levels below a given wavenumber k or 
energy E, in a desymmetrized billiard, is given by the 
Weyl formula N = k 2 /16 = E/8. 

The effective spread in angular momentum of an eigen- 
states \& is characterized by the localization length I 
which measures the typical number of p n (l) values which 
are substantially different from zero, or the size of the 
angular momentum interval in which p(l) is significantly 
different than zero. Of course, there is no a unique def- 
inition of localization length. In particular the choice is 
quite delicate for the quantum stadium where localiza- 
tion is algebraic unlike billiards with smooth boundary 



where localization is expected to be exponential. We have 
therefore considered four different working definitions: (i) 
spread in angular momentum 

£ 1 =/3 1 ((v[/|f 2 |*)-(vl/|f|vl;)2)i/2 j 

(ii) information entropy 

£ 2 = /3 2 exp^-5>(Z)lnK0) 

(iii) inverse participation ratio, 

i 

(iv) 99% probability localization length 

£ A =p 4 imn{#A;J2pk(l) > 0.99}. (16) 

The numerical constants f3j, j = 1,2,3,4 have been deter- 
mined with the condition that £j should take the max- 
imal value Z max = k — \/2E in the regime of quantum 
crgodicity. All the four definitions of localization length 
give qualitatively the same results. However, the 99% 
localization length (|l6| ) which is proportional to the min- 
imal number of angular momentum eigenvalues that are 
needed to support 99% probability, is the least sensitive 
to the slowly (algebraically) decaying tails of the distribu- 
tion p(l) and has given the sharpest and less fluctuating 
numerical results. The results shown in figures 3 and 
4 were obtained using the definition ([r]) of localization 
length with the numerical value/34 = 1.38. 

C. Computation of Husimi functions 

In figure 5 we show Husimi functions of (quarter) sta- 
dium eigenstates on the Poincare-Birhkoff section, with 
coordinates (<p,l)- The desymmetrized phase space is a 
rectangle ((p,l) G [0,7r/2] x [0,k],k = V2E. We consider 
(almost) minimal wavepackets \W) at phase space point 
(8, 1), which in the angle representation \ip), < (p < 2tt 
read 

(tp\0,l) = Cexp(-^-sm 2 (^-6))+il^\ (17) 

where C is a normalization constant and a is the aspect 
ratio. We take a = 1 which means that the wavepacket 
has the same relative size in both directions. The Husimi 
function h(ip, I) of a billiard eigenstate \&(r) is computed 
w.r.t. the normal derivative of an eigenfunction along the 
boundary r(ip), d„^ = (df/dip) A W/\\dr/dtp\\, 

h(9,l) = | J d(t>^\0,l)d n ^(f^))\ 2 , (18) 

where the numerical integration along the boundary of 
the billiard has been performed accurately using high- 
order Gaussian quadratures. 



VIII. CONCLUSIONS 

The billiard in a stadium is a very clean mathematical 
model of classical chaos and it is therefore particularly 
convenient to study the modifications that quantum me- 
chanics introduces in our general picture of deterministic 
chaos. While the classical motion is ergodic and mixing 
for any non zero values of the parameter e, the quantum 
motion, as we have seen in this paper, exhibits a very rich 
structure and different regimes of motion as a function of 
the parameter e or the energy E. Two are the main points 
we would like to stress: a) classical cantori may have 
strong effects on quantum mechanics. They lead to a new 
quantum border which is distinct from the perturbative 
border and from the localization border. In the regime 
of quantum cantori the rescaled localization length tjk 
does not depend on energy or wavenumber k. In other 
words, the quantum dynamics is basically determined by 
the classical structure. This effect should be observable 
in real experiments, b) The mechanism of localization is 
strictly connected to the sparsity of EF's when expanded 
on the basis of unperturbed circle states (or vice versa, 
sparsity of unperturbed states when expanded in terms 
of the stadium EF's). When we increase perturbation 
(or decrease h), sparsity decreases, localization length of 
EF's increases until the delocalized or quantum ergodic 
regime is reached. We conjecture that this mechanism 
of localization is typical of conservative Hamiltonian sys- 
tems. Quite interestingly, the mechanism of localization 
in WBRM is quite different and is associated to an asym- 
metry between EF's (solid and narrow) and unperturbed 
states (delocalized but sparse). In this sense, WBRM 
cannot be taken to represent the typical behaviour of 
classically chaotic conservative systems. 

The main impulse to the present paper originated from 
a previous work of B.Chirikov |25| in which a wiggled cir- 
cle billiard was introduced as a model for localization in 
conservative systems. Also we have benefit from several 
results on the standard map. This confirms the impor- 
tance and generality of the Chirikov standard map for the 
study of classical and quantum deterministic chaotic mo- 
tion. Since 30 years, generations of physicists all around 
the world, have learned the beauty and complexity of 
nonlinear dynamics on the base of this apparently simple 
map. We believe that the Chirikov map is at the root 
of the impressive growth of the whole field of nonlinear 
dynamics and chaos. Strangely enough this fundamental 
contribution of B.V. Chirikov to one of the main achieve- 
ments of physics of this century has never been properly 
recognized. Several years ago, in a friendly private letter 
to a colleague concerning a paper in which he failed to 
give the necessary credit to Chirikov, a common friend 
Joseph Ford, now deceased, wrote " ..why not give to the 
old Russian bear his due?". We would like to turn this 
question to the scientific community. One of the authors 
of the present paper(GC), would like to take the opportu- 
nity of this special issue of the journal to express his deep 
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FIG. 1. The rescaled classical diffusion rate (4) 

D(l)/{e 5/2 (2E-l 2 )} asafunction of the local value L = l/y/2E 
of the angular momentum for different values of the control 
parameter e. 




0.16 
L=l/k 

FIG. 2. Time evolution of a single classical orbit, followed 
up to 20, 000 bounces, for the classical billiard with e = 0.003. 
The orbit is initially started in the middle of the largest 'is- 
land' ( L = 0,9 = tt/4- 0.0016) (fig 2a). Angular momentum 
probability distribution p(l) of the corresponding eigenstate 
with e = 0.003 and eigenvalue k = 5999.8166. As it is seen the 
state is uniformly distributed over the cantorus in the main 
island (Fig. 2b). Notice also the same eigenstate in Husimi 
phase-space representation (2nd state in Fig. 5b). 
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FIG. 3. Rescaled localization length a — l/k versus the 
scaling variable x = e 3 ^ 2 k for five different values of e 
(60 < k < 12, 000). Each point is obtained by averaging over 
a large number v of consecutive eigenstates (y = 100 for small 
k and v = 1, 000 for large k). The numerical data clearly show 
the cantori border x — 1. In the cantori region l/k is con- 
stant as expected, while for x > 1 the numerical data agree 
with the theoretical prediction ©(dotted curves). For large 
x, the value of localization length a approaches the maximal 
ergodic value a = 1. 
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FIG. 4. Rescaled localization length a versus e for two dif- 
ferent wavenumbers k. The full lines give the classical esti- 
mate for the average width of cantori. It is seen that be- 
low the cantori border x = l(e = 0.003 for k — 6,000 and 
e — 0.0063 for k — 2, 000) the localization length is propor- 
tional to e and independent on k. Above the border x = 1 
instead, the numerical data follow the analytical estimate (15) 
(dotted curves). 



Color figure 5 (a-d) in (Windows/OS2 bitmap format) 
can be obtained by e-mail request on prosen@fiz.uni-lj.si 



FIG. 5. The Husimi phase space representation of quan- 
tum eigenstates of the stadium for different values of pa- 
rameters e and k — y/2E. The Husimi function of a 
normal derivative of an eigenfunction along the boundary 
of desymmetrized (quarter) stadium billiard (section VII) 
is plotted (x-axes: angle < <p < ir/2, y-axes: angu- 
lar momentum < I < \Pl~E). We show four consec- 
utive eigenstates in the following order: up-left, up-right, 
down-left, down-right for (a) e = 0.001, k = 5999.9775, 
5999.9812, 5999.9815, 5999.9818, (N « 2.25 • 10 6 ), (b) 
e = 0.003, k = 5999.8138, 5999.8166, 5999.8165, 5999.8175, 
(c) e = 0.01, k = 6000.0079, 6000.0104, 6000.0112, 6000.0121, 
and (d) e = 0.01, k = 12000.0030, 12000.0040, 12000.0041, 
12000.0051, (N « 9 • 10 6 ). Note that the second state in (b) 
is the same as in figure 2b. The color scale (to the right of 
each figure) is proportional to the logarithm of Husimi func- 
tion (white is large, black is small, entire color scale extends 
over 8 orders of magnitude, a factor of 10 . The cases (a,b) 
are in cantori region and one should note that the sizes of 
bright regions are roughly the same, while (c,d) are in the 
region of dynamical localization. 
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FIG. 6. Average shape of eigenstates in angular momen- 
tum representation, (before averaging all the eigenstates have 
been shifted to the same mean value (I) of angular momentum. 
The average has been made over about thousand consecutive 
eigenstates for e = 0.005, k £s 12000 (regime of dynamical lo- 
calization), and e = 0.001, k « 6000 (cantori regime). The 
dotted lines give the expected slope —4 for the tails. 
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FIG. 7. The nearest neighbour energy level spacing 
distribution P(S). We show data for the two spectral 
stretches in the same energy window, 1993.8 < k < 2006.2 
(N w 250,000), containing about 3,100 levels each, namely 
e — 0.002 (cantori region), and e = 0.02 (region of dynamical 
localization). The theoretical Poisson and Wigner statistics 
are shown for comparison. Notice that in the cantori region 
e = 0.002 (much above the perturbative border) the numeri- 
cal data closely follow the Poissonian statistics (the deviation 
is only slightly larger than the expected statistical error) . 
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FIG. 8. The structure of the matrix of coefficients 
c™ m in the regime of dynamical localization: e = 0.01, 
k ~ fco := 2000.0 (Fig.8a) A small dot at abscissa k° am — ko 
and ordinate k n — fco is plotted if |c™ m | > 0.02 and a large 
dot if |c™ m | 2 > 0.1. Notice that the total number of levels 
along each axis is about 1250. In Fig. 8b we plot the average 
local density of states w(k — ko) and the average eigenfunction 
W (k — ko) for the same e = 0.01, averaged over a stretch of 
5000 consecutive eigenstates in the interval 1990 < k < 2010. 
The dotted curve gives the theoretical LDOS for an ergodic 
billiard W e (k) = (sin(fc)/fc) 2 /7r. 



